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Abstract 

Uncertainties in the next-to-leading-order calculations of heavy- 
quark {Q) production are investigated. Predictions for total cross 
sections, single-inclusive distributions of heavy quarks and heavy 
flavoured hadrons, as well as for QQ correlations, are compared 
with charm and bottom data. 

The description of heavy-quarkonium production requires a sep- 
aration of the short- distance scale of QQ production, which is set 
by the heavy-quark mass from the longer- distance scales associated 
with the bound-state formation. Various factorization approaches 
are compared, in particular with respect to the different constraints 
imposed on the colour and angular-momentum states of the QQ 
pair(s) within a specific quarkonium state. Theoretical predic- 
tions are confronted with data on heavy-quarkonium production 
at fixed-target experiments and also at pp colliders, where frag- 
mentation gives the leading- twist cross section in l/p\ and l/m? 
at high transverse momentum. 
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1 Production of heavy quarks 



1.1 Theoretical status 

Total cross sections and single-inclusive distributions of open heavy-quark pro- 
duction have been calculated in perturbative QCD (pQCD) to next-to-leading- 
order (NLO) accuracy for essentially all high-energy reactions, hadroproduction 

0. photoproduction @], leptoproduction (deep- inelastic lepton-nucleon scatter- 
ing) 1^, 77 coUisions [Q, and deep-inelastic 67 scattering In the case of 
photo- and hadroproduction, the NLO calculations have even been implemented 
in a Monte Carlo programme such that also double-differential distributions 
can be studied. Attempts are ongoing to develop such an exclusive programme 
also for deep- inelastic scattering 0. 

In pQCD, the differential cross section of heavy-quark (Q) production in the 
collision of hadrons A and B is given by the factorized expression 

da[AB ^QQX]{pa,Pb) = 

J2 J dxi dx2fi/A{xi, fJ^p) fj/Bix2, f4,)da[ij -> QQX']{xiPa, X2PB, fJ^F, IJ^r) ■ (1) 

id 

Here i, j represent the interacting partons (gluons, light quarks and antiquarks) 
and the functions fi/A{x,^p) are their number densities, the parton distribu- 
tion functions (PDF) evaluated at momentum fraction x and factorization scale 
fip- The short-distance cross section a is calculable as a perturbation series in 
c^sif^R) where the strong coupling constant is evaluated at the renormalization 
scale fiR- Leading-order (LO) diagrams for heavy-quark hadroproduction are 
shown in Fig. 1. 

The obvious question is, of course, whether pQCD does describe the exper- 
imental data on charm and bottom production. To answer this question the 
uncertainties in the theoretical predictions have to be investigated. These arise 
from four main sources: 

1. The heavy-quark mass rriQ. Conservative ranges are 1.2 < m,c < 1.8 GeV and 
4.5 <mb< 5.0 GeV. 

2. The size of unknown higher-order corrections. An estimate may be obtained 
by varying the factorization scale fip and the renormalization scale fiR around 
their 'natural' value fiR = fJ^F = say between m/2 and 2m. In the case of 
the charm quark, one has to note that most PDF /i/H(x,/i|n) are not valid for 
/ii? < 2 GeV and that the running of «s(/i|j) is no longer purely perturbative for 

< mc/2. Proper error estimates are therefore not as easy as in the case of the 
h quark (see below). 

3. The value of the QCD scale A and the shape of the PDF. Since these are 
strongly correlated (in particular A and the gluon density), the choice of differ- 
ent A values (in 60 < A5 < 300 MeV, say) in the partonic cross sections should 



1 



ideally be accompanied by parametrizations of PDF that have been fitted with 
those values of A. However, LEP data (and also the r hadronic width) favour A 
values that are almost a factor of two larger than is found in fits to deep-inelastic 
scattering data. In turn, one does not cover the full range of uncertainty in A 
when restricting to the available PDF parametrizations. In particular, the cross 
sections may be larger since larger values of A imply larger values of Ug- In the 
absence of fits with A frozen to the desired values one may choose to simply 
change the value of A in the partonic cross sections. 

4. Non-perturbative effects such as heavy flavoured hadron formation and the 
intrinsic kj- of the incoming partons. 



1.2 Low-energy total cross sections 

In the case of charm production it turns out that the biggest uncertainty arises 
from the variation of m^. In ref. it is found that predictions of the lowest 
and largest values of the total charm cross section in hadroproduction may differ 
by up to a factor of 10 from the central value (i.e. the highest value is a factor 
of 100 larger than the lowest one)! Predictability substantially improves when 
increasing the quark mass: for h quarks, the theoretical error is reduced to a 
factor of ±2-3. Both for charm and for bottom, the predicted total cross sections 
are in reasonable agreement with recent data (Fig. 2); the charm data are claimed 
to be consistent with a charm-quark mass of 1.5 GeV f^. 

The uncertainty bands of the theoretical predictions in Fig. 2 represent the 
maximum variation of the cross sections in a parameter range defined as follows 
[||. The scale hr was varied in [mQ/2,2mQ] and A4 in [100,300]MeV in order 
to account for a range of uncertainty in A wider than is provided by PDF fits. 
However, only two particular sets of PDF parametrizations (one for the proton 
and one for the pion 0) were chosen with A4 = 190 MeV. Hence the correlation 
between A and the (nucleon and pion) PDF was not properly taken into account. 
The mere variation of A overestimates the associated uncertainty. 

The scale was also varied in [m;,/2,2mf,] in the case of bottom, but kept 
constant at 2mc in the case of charm since the chosen PDF parametrizations 
are not valid at smaller scales. The choice = maximizes the values of the 
PDF for most values of x and in turn of the charm cross sections. Moreover, PDF 
decrease (eventually non-perturbatively) as — >• 0. Finally, the lowest allowed 
value of iiR (m™"/2 = 0.6 GeV), together with the largest value of A (300 MeV), 
gives an as value of 1.21, which is certainly beyond QCD perturbation theory. 
In conclusion, the error bands in ref. are fairly pessimistic, in particular for 
charm, and values of smaller than 1.5 GeV are favoured by the data. 

The question whether the NLO calculation can accommodate the hadropro- 



duction data for rric = 1.5 GeV was also addressed in ref. |T0[. The cross section 



for the choice hr = hf = rric = 1-5 GeV (and the PDF parametrization GRV |]TT 
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that extends down to fip ~ 0.6 GeV) lies a factor of 2-3 below the data (Fig. 3). 
A word of caveat on the data shown in Fig. 3 on the total cc production cross 
section from pp and pA interactions is in order. 

First, in fixed-target experiments mostly single charmed mesons in the region 
xp > are measured. Additional, poorly known information is needed to con- 
vert a{D/ D)[xp > 0] into a[cc]. So far only premature extractions of a[cc] exist, 
obtained by simple assumptions on cr(Ac), cr{Ds), the xp distributions outside 
the measured range, in particular in pA collisions, the distortion of the xp distri- 
butions due to non-perturbative effects, etc. This explains (at least partly) the 
spread in the low-energy {^/s < 40 GeV) data shown in Fig. 3. 

Second, the ISR data (53 < ^/s < 62 GeV), inferred from lepton measure- 
ments in coincidence with a reconstructed charmed hadron, are presumably too 
large due to the assumed shape of the production cross sections: flat distribu- 
tions in Xp for the A^ and (1 — xp)^ for the D. Restricting to the most recent 
experiments only consistency is found P, |1^ between the data (Figs. 2 and 3) 



and these data agree rather well with fiR = fip = rric = 1.3 GeV in the NLO 



calculation using the GRV PDF ||TT|. The result is similar (Fig. 3) with the MRS 
distributions [l^ and fiR = fip = 2mc where now the large fiR requires an even 
smaller charm-quark mass, rric ~ 1.2 GeV. Note that such small values of rric sug- 
gest that the bulk of the total cross section comes from invariant masses smaller 
than 2m D- 

The theoretical uncertainty is smaller in photoproduction: there the Born 
cross section starts at 0{as), compared to O(a^) in hadroproduction, and also 
only a single parton density enters. For example, the charm cross section can now 
be calculated up to a factor of about 3 |^ . The theoretically better-controlled b- 
quark cross section has unfortunately not yet been measured in photoproduction; 
the forthcoming measurement in 'jp collisions at HERA is eagerly awaited. 

Additional tests of heavy-quark production will soon also be possible in elec- 
troproduction of heavy quarks at HERA and in electron-photon collisions at 
LEP2. The most reliable predictions can be made for heavy-quark production in 
77 collisions. Even for c quarks, the total uncertainty is only about ±(30-35)% 
[0. The experimental situation is so far controversial: taken at face value, the 
cross section seems to fall when going from Tristan |[l^ (v^ee ~ ^0 GeV) to LEPl 
H (v^,e~90GeV)! 



1.3 Single-inclusive distributions 

Single-inclusive distributions can be calculated in pQCD by folding the (NLO) 
cross sections of heavy-quark production with the corresponding (NLO) frag- 
menting function, e.g. c ^ D + X . In general, the xp and pp spectra of the 
heavy-flavour hadron will be softer than those of the mother heavy-quark. In the 
case of charm hadroproduction, the measured D-meson pp distribution is con- 
siderably harder than the one calculated in NLO perturbation theory. Since the 
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c.m.s. energies are small (^i = 20-40 GeV), it appears unlikely that the discrep- 
ancy can solely be resolved by yet higher-order corrections. Rather the description 
of differential charm distributions requires the inclusion of non-perturbative ef- 
fects. These may arise from the intrinsic transverse momenta kx of the incident 
partons and/or higher- twist corrections in the hadronization. 

The effect of intrinsic kr is well established in the Drell-Yan process and 
its inclusion gives also reasonable agreement |^ for D-meson distributions. 
In fact, photoproduction data are reproduced with kx values as they are also 
observed in the Drell-Yan process, i.e. {k^Y^"^ ~ 0.7 GeV. On the other hand, 
hadroproduction data require values that are as large {{k^y^'^ ~ 1.4 GeV) as the 
charm-quark mass. Further studies are certainly needed, in particular on the 
strong correlation of kx with (larger yields harder px spectra and, in turn, 
allows smaller fc^ values). 

Distributions in xp have so far only be measured for charmed hadrons and 
do not at all agree with the pQCD predictions (they are too soft), even after 
inclusion of the kx effects . Indeed, an agreement must not be expected in the 
first place, as factorization, underlying the perturbative approach, is valid only 
at large pt and hence breaks down at large xp. Additional non-perturbative con- 
tributions associated with the hadronization become important: The enhanced 
production of D-mesons whose light valence quark is of the same flavour as one of 
the valence quarks in the beam hadron (leading-particle effect) and the dragging 
of charm quarks in the colour field of the beam fragments (colour-drag effect). In- 
deed, a fragmentation model that contains these effects, such as the Lund string 
fragmentation one, successfully describes the observed spectra. An alternative 
explanation |]T3|] invokes new production mechanisms for hard processes at large 
xp through terms that are formally higher-twist effects oc I/ttlq but are enhanced 
by inverse powers of 1 — xp- 

1.4 QQ correlations 

In leading order, pQCD predicts that the heavy quark and heavy antiquark are 
produced exactly back-to-back (Fig. 1), implying priQQ) = and A0 = vr, 
where Pt{QQ) is the transverse momentum of the pair and A0 the angle between 
the projections of the momenta of the pair onto the transverse (w.r.t. beam 
axis) plane. NLO corrections, as well as non-perturbative effects, can cause a 
broadening of these distributions. The question whether NLO predictions can 
account for the available experimental data in hadro- and photoproduction of 
charmed particles was addressed in ref. f^. Agreement with data was found for 
the A0 distribution after inclusion of a modest intrinsic transverse momentum 
(fc2)V2 _ 0.7-1 GeV. 

The theoretical pxiQQ) distributions, on the other hand, are too soft unless 
much larger kp values are used. Although it is mostly the intrinsic kp of gluons 
that enters heavy-quark production and not the kp of quarks that is known from 
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the Drell-Yan process to be of the order of 700 MeV, it yet appears unreasonable 
to have a that is of the order of the hard scale, namely ~ mg. Three possible 
explanations come up. First higher-order corrections can broaden the Pt{QQ) 
shape. After all, the O(a^) calculation, although being of NLO accuracy for the 
total cross section and single-inclusive distributions, is still of leading order for 
the A(f) and priQQ) ones. Second, further non-perturbative contributions besides 
kx could be significant. And third, experimental high-pj^ (cc) data are still rather 
limited in statistics and might possibly come down. New measurements of charm 
hadroproduction by WA92 are on their way . The different hypotheses to explain 
the high-pT^QQ) discrepancy can also be tested in photoproduction at HERA 
where high-px^QQ) data should soon become available. 



1.5 Bottom production at collider energies 

Owing to its heavier mass, not only the total b cross section but also differential 
distributions of the 6-quark and 6-fiavoured hadrons can be predicted with higher 
accuracy (at least, if mg/y^ is not too small). In fact, the pt distribution 
measured at both 630 GeV and 1.8 TeV is now found to be consistent (Fig. 4) 
with the NLO predictions 0]: On the one hand, the experimental data from the 
Tevatron 'came down' through improved b tagging, 6 — > J/tp, ip' via secondary b 
vertices, B J /if) + K^*\ bb correlations, b — > fxX identification in DO. On the 
other hand, the theoretical prediction 'came up', partly because PDF are now 
known to steeply increase in the relevant x range |jTB| . 

The biggest change in ref. |^ compared to previous estimates is, however, 
that A is allowed to be as large as the central LEP value (A5 ~ 300 MeV or 
A4 ~ 420 MeV). The upper curve in Fig. 4 is obtained by stretching all parameters 
to their extremes: small quark mass [nib = 4.5 GeV), small renormalization scale 

il^R = I^F = \Jpt + ^6/2), and large QCD scale (A5 = 300 MeV). Yet, at the 
same time the MRSA PDF parametrization [|1^ with A5 = 151 MeV is being 
used. First estimates [Q indicate that varying A within a limited range, without 
refitting the PDF, is not as large an overestimate of the systematic effect of the 
A uncertainty on the b cross section as one might expect at first. Given the 
potential sensitivity of the b cross section on small-x effects, further theoretical 
studies should certainly be pursued. 



2 Bound-state production 

The production of quarkonium states below the open charm/bottom thresholds 
presents a particular challenge to QCD. Because of the relatively large quark 
masses, c and b production are perturbatively calculable, see previous section. 
However, the subsequent transition from the predominantly colour-octet QQ pairs 
to physical quarkonium states introduces non-perturbative aspects. Three models 
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of bound-state formation are discussed in the literature. In order of increasing 



sophistication, these are the colour-evaporation model |T^, |T9l , the colour-singlet 



model ||2y, |T9|, and the one based on a new factorization formula. 

In QCD, the total angular momentum J, the parity P, and the charge conju- 
gation C are exactly conserved quantum numbers. Hence the energy eigenstates 
\H) of heavy quarkonium can be labelled by the quantum numbers J^*-^ (besides 
a quantum number n = 1, 2, . . ., to distinguish states identical apart from their 
mass M as the 1 states J/ip and ip'). Phenomenologically very successful is 
the non-relativistic potential model, in which \H) is assumed to be a pure quark- 
antiquark state \QQ)- Obviously, the QQ pair must then be in a colour-singlet 
state 1 and in an angular-momentum state "^^^^Lj that is consistent with the 
quantum numbers J^'" of the meson (P = (—1)'^+^, C = (—1)'^^'^) 

potential model : \H{nJ^^)) = IQQin'^^^^Lj,!)) . (2) 

Here S* = 0, 1 is the total spin of the quark and antiquark, L = 0, 1, 2, . . . or 
S, P, D, . . . is the orbital angular momentum, J is the total angular momentum, 
and n denotes the principal quantum number. 



2.1 Colour-singlet model 

In the colour-singlet model, the dominant production mechanism of a heavy 
quarkonium is assumed to be the one in which the quarkonium is produced at 
short distances in a colour-singlet QQ state with the correct quantum numbers. 
Hence the cross section is given by the factorized form {v = xp, pt, ■■■)'■ 

da[H{nr^)]{s,v) ^ ^ ^ da[QQin''+'Lj,l)]{s,v) .3. 
d'u " dv 

The non-perturbative probability Fni for the QQ pair to form the bound state H 
is given in a calculable way in terms of the radial wave function or its derivatives 



IP 



nL 



FnL (X ' (4) 



M: 



3+2L 
H 



and can either be calculated using a phenomenological potential or extracted 
from the H decay widths that are given by an expression similar to ||19|| . 

Actually, (|]) gives the dominant cross section only if the relevant momen- 
tum scale Q is of the order of the heavy quark mass. An example is the total 
hadroproduction cross section of J/ip where ~ s ~ 4m^. At large scale Q, 
however, the short- distance production becomes suppressed by a factor im? /Q"^ 
with respect to production via fragmentation that then gives the leading-twist 
cross section in l/Q"^ and 1/m?. A long-known example ETI] is J/ip production 



in Z° decays where cc pair production followed by the fragmentation c ^ J/ip 
dominates over J/ipgg. 
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Similarly, fragmentation processes start to dominate pp H + X at high pt 
. In this case, important contributions arise from gluon fragmentation 

da[pp H{p) + X]{s,pt) _ 
dpT 

dz Dg^H[z,ii,m) . (5) 

dpT 

The fragmentation functions Da^niz, fJ'jfn) specify the probability for partons a 
(gluons, light and heavy quarks) to hadronize into the hadron if as a function of 
its longitudinal momentum fraction z relative to a. Large logarithms of pt//^ in 
the parton cross sections a are avoided by choosing the (arbitrary) factorization 
scale n of the order of the large scale pt- Large logarithms of fi/m then necessarily 
appear in the fragmentation functions Da^niz, fi,m), but they can be summed 
up by evolving the "input" distributions Df\^{z) = Da^H{z,m,m) with the 
standard evolution equations. 

In the colour-singlet model, the fragmentation functions at the input scale 
m, Df}^^{z) can be calculated as a series in ^^(m) by assuming that they take 
the same factorized form as (|^). For example, the lowest-order diagrams that 
contribute to gluon fragmentation into J/ip are g cc(^Si,l)gg, so that 



Dflj/A') = (^]\Rism'm + 0{at) , (6) 



where f{z) is a calculable function. Note that the kinematic regime in which 
the a ^ H fragmentation graphs become important occurs when the lab-frame 
energy Ea of the parton a is large, but its squared four-momentum is close to 
the square of the charmonium bound-state's mass ~ 4m^. Terms subdominant 
in the ratio p^/E"^ are therefore neglected. 

In contrast to 5- wave fragmentation functions |2^, |2^, the colour-singlet 



contributions to fragmentation functions into xqj state are, however, singu- 
lar. The process g — >■ cc{^Pj,l) +g diverges logarithmically (analogously 



;i - z)+ eo 



777 

+ \n-5{l-z) 



+ rjiz) 



c^cc(3pj,i) + cii) 

'(7) 

Here the infrared divergence, associated with the soft limit of the final-state gluon 
has been made explicit by the introduction of a lower cutoff eo on the gluon energy 
in the quarkonium rest frame. The presence of the infrared sensitive term clearly 
spoils the factorization assumption of the colour- singlet model. 

In order to still separate the long- and short-distance contributions (at least) 
a second non-perturbative parameter has to be introduced 

Dfl^^iz) = ^^^^0, + dsiz) 08(A) . (8) 
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The first term in describes the emission of a perturbative gluon with energy 
above some cutoff A 



81 m-^ 



(1-^)+ A 



ln^5(l-z) 



+ rj{z)\ .(9) 



In the second term 



16 A (2J+l)gi 
OslA) = ^^a,ln- — 



27t[ eo 
24 



C^8(^) = TT7 ^^'^(l-'^) ' (10) 



the presence of the infrared scale eo indicates that Os actually has to be considered 
as an additional non-perturbative parameter besides R'p{0), as expression (p!0|) 
for Os can at best give the perturbative part of Os- 

Even with improved P-wave factorization, the colour-singlet model still fails 
in the description of data on high-p^ prompt charmonium production. Such data 
now become available at the Tevatron where vertex detectors can be used to 
separate the charmonium states coming from b quarks from those that are pro- 
duced by QCD interactions. Through the inclusion of fragmentation mechanisms 
the theoretical predictions for prompt J/ip production can be brought to within 
a factor of 3 of the data ^Q, 27], close enough that the remaining discrepancy 



might be attributed to theoretical uncertainties (Fig. 5). 

However, this conclusion of a successful description of J/ip production relies 
on the postulation of a very large g Xcj{^P) fragmentation contribution^ where 
the Xcj subsequently decay into ■yJ/ip- Since such a contribution is absent for ip' 
the prediction of its production rate falls a factor of 30 below the data (Fig. 5), 
casting doubtsQ on whether the J/ ip production is correct at all. These doubts 
are strengthened by the observation that the ratio of ip' to J/ip measured at 
high transverse momenta at the Tevatron is quite compatible ^ with the 
Pt- integrated fixed-target and ISR data where direct J/ip production is known 
[ P^ to dominate the indirect production via Xcj decays. Hence the major part 
of prompt high-pj- J/tp^s should be directly produced rather than originate from 
radiative Xcj decays. 

The breakdown of the simple factorization into a single long-distance matrix 
element and a short- distance Wilson coefficient, discussed above for gluon frag- 
mentation into P-wave quarkonia, has, in fact, previously been pointed at for 

^ This large contribution follows from an optimistic choice of Os in (§)■ 
^The ad hoc introduction of (not-yet) observed charmonium states above the DD threshold 
such as higher P-wave or Z?-wave states appears quite questionable. Only with very optimistic 
production rates and branching fractions into ip' can one account for the observed rate for 
prompt ip' production pq]. 
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the cases of _B-meson decays into Xcj ^T[l' hadronic xqj decays, and total cross 
sections of xqj hadroproduction ||l9l. Also these analyses showed that, even if 



the colour-singlet model is extended through the inclusion of colour-octet mech- 
anisms for P-wave states (analogues of the Os term in (||), see also section |2.3|) , 
the model cannot accommodate for all the data on quarkonium production and 
decays. 

It must therefore be concluded that also in the case of S'-waves, a generaliza- 
tion of the naive factorization a la (0,^) is necessary, at least in cases where the 
colour-singlet mechanism is suppressed by a short- distance coefficient, i.e. by a 
high power of asijn), cf. (^. A factorization of long- and short-distance physics 
more generalized than the one assumed in the colour- singlet model will be pre- 
sented in section p.3| . First, however, a quarkonium production-model diametrally 
opposite to the colour- singlet model is discussed. 



2.2 Colour-evaporation model 

The starting, i.e. perturbatively calculated cross section for heavy quarkonium 
production in the colour-evaporation model is the (usual) cross section of open 
heavy-quark pair production that is summed over all spin and colour states. 
Hence all the information on the non-perturbative transition of the QQ pair to 
the heavy quarkonium H of quantum numbers J^^ is contained in "fudge factors" 
F^jpc that a priori may depend on all quantum numbers 

da\H{njP^)]{s,v) dd[QQ]{s,v) 

Tv = Tv • ^^^^ 

In (0) (y[QQ\ is the total "hidden" cross section of (open) heavy-quark produc- 
tion calculated by integrating over the cc pair mass (in the case of charmonium) 
from 2mc to 2m£,. For example, in hadronic collisions at high energy, the domi- 
nant production mechanism is gluon fusion, so that (cf. (|I])): 




fglA{xi,m^)fglB{x2,m^)a[gg^ccX]{s)8{s-xiX2s) . (12) 

Note that a[cc\ is the spin-summed cross section and that the heavy-quark pair 
can be both in a colour- singlet and a colour-octet state (in LO, qq annihilation 
produces only colour-octet cc pairs, while gg fusion also leads to colour-singlet 
states). The colour-octet cross section is the dominant one. The cc configuration 
arranges itself into a definite outgoing charmonium state by interacting with the 
collision- induced colour field ("colour evaporation"). During this process, the c 
and the c either combine with light quarks to produce charmed mesons, or they 
bind with each other to form a charmonium state. 
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As shown in ref. more than half of the subthreshold cross section cr[cc] in 
fact goes into open charm production (assuming < m[ri(\/2 ~ 1.5 GeV); the 
additional energy needed to produce charmed hadrons is obtained (in general 
non-perturbatively) from the colour field in the interaction region. The yield of 
all charmonium states below the DD threshold is thus only a part of the total 



sub-threshold cross section: in this aspect the modern version of the model |^ is 



a generalization of the original colour-evaporation model ||18|, pi)[|, which neglected 
the contribution of a[cc] to open charm production. Using duality arguments, it 
equated a[cc\ to the sum over the charmonium states below the DD threshold. 

Neither the division of a"[cc] into open charm and charmonia nor the relative 
charmonium production rates are specified by the generalized colour-evaporation 
model. Hence its essential prediction is that the dynamics of charmonium pro- 
duction is that of a"[cc], i.e. the energy dependence, xp- and p^-distributions of 
H are identical to those of the free cc pair. In particular, ratios of different char- 
monium production cross sections should be energy-, xp-, and pT-independent. 
In other words, the non-perturbative factors F^jpc should be universal constants 
whose values may, however, depend on the heavy-quark mass. In contrast to 
earlier expectations [0, a recent comprehensive comparison of the generalized 
colour-evaporation model with data indeed confirmed these expectations ||30| . 

Figure 6 shows the ratio of J/V' production from the decay Xc ^ 1 J/'ip 
to the total J/ip production rate, which provides a measure of the Xc/{J/4') 
rate, and Fig. 7 shows the measured ifj' / {J / ip) ratio. Both ratios are found to 
be independent of the incident energy, the projectile (pion or proton), and the 
target (from protons to the heaviest nuclei). Moreover, it is noteworthy [^, ^ 



that the ratio if)' /{J/ if)) measured at high transverse momenta at the Tevatron is 
quite compatible with the pr-integrated fixed-target and ISR data (Fig. 8). Also 
the available bottomonium data agree with constant production ratios ||30| . 

Figures 9 and 10 show the energy dependence of J/ip and T production in 
hadronic collisions; the agreement with data over a wide range is rather impres- 
sive. Because the data generally give the sum of T, T' and T" production, the 
measured cross section for the sum of the three T states in the dilepton decay 
channel is shown, denoted by B{da/dy)y=Q. Using the average values of the T" to 
T and T' to T production ratios, 0.53 ± 0.13 and 0.17 ± 0.06, respectively, fits 
to the data in Figs. 9 and 10 yield 



2.5 X 10 ^ charm 

4.6 X 10-2 bottom 



'11 ~ /I « XX in-2 K^^-+^v^ ■ (-'-^) 



The hidden heavy- flavour cross sections &[QQ] in Figs. 9 and 10 were calcu- 
lated in NLO using the MRS D— ' parametrization [|T2 of PDF with fir = fJ'R = 



2mc = 2.4 GeV and fxp = jjR = rrib = 4.75 GeV, respectively. These parame- 
ters provide an adequate description of open heavy- flavour production [|l^, cf. 



section 1. Results similar to (|1^) are obtained if one uses other choices of the 
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parameters that are tuned to the open heavy-flavour data, for instance the GRV 

1.3 GeV, and 



Tl[] with fip = fifi = niQ, 



rrir 



4.75 GeV. 



parametrization 

From ([T3| ) one concludes that the fraction of a[cc] producing charmonium rather 
than open charm is about 10%. 

Equally good agreement is found for the energy dependence of J/ ip production 
with pion beams. However, the fraction of J/ ip in the hidden charm cross section 
must be slightly higher to reproduce the pion data well, with Fn — = 0.034 for 
a good fit. This may well be due to greater uncertainties in the pionic parton 
distribution functions. Consistency with the colour-evaporation model is also 
found for the longitudinal momentum dependence of charmonium production. 
The calculations for the xp dependence of J/ip production agree with data 
(Fig. 11) from low-energy pp interactions (where qq annihilation is important) up 
to pp interactions at the highest available energy and out to the largest xp values 
probed (where the "intrinsic charm" mechanism could have been important). 

Concerning the pt distribution, the model provides essentially no prediction 
for low-px charmonium production. There is the intrinsic transverse momentum 
of the initial partons, the intrinsic momentum fluctuations of the colour field 
which neutralizes the colour of the cc system in the evaporation process and, 
at larger pt, higher-order perturbative terms. Since there is no way to separate 
these different contributions in the low-pT region, the model has no predictive 
power. 

On the other hand, the high-p^ tail should be successfully describable in the 
colour-evaporation model, presumably with the same normalization (0). The 
time scale to form a quarkonium bound state is much larger than the one to 
produce the (compact) cc pair. Hence the fraction of the cc cross section (in 
2mc < Mcc < 2m£)) that becomes a. J/ip (or ip') should be independent of the 
cc production, i.e. the same for low-pT and high-p^ processes. By the same 
argument, the universal ratio of ip' to J /ip production, observed experimentally, 
can depend only on the relative magnitude of the respective wave functions at 
the origin 



e ' e 



M 



am 



r(W 



(14) 



a{J/ij) 



exp 



Relation (|T^) holds to very good approximation p9| , pO[| . Note, finally, that gluon 
fragmentation gg g*g with g* ^ cc ^ J/ip is part of the lowest-order (0(a^)) 
diagrams describing high-p^ charmonium production in hadronic collisions, while 
charm fragmentation gg c*c with c* ^ J/ip first occurs at O(a^). 
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2.3 A new factorization approach 

The two previous sections discussed two extreme scenarios to describe production 
cross sections (and decay rates) of heavy quarkonia. In the colour-evaporation 
model (p!TD, no constraints are imposed on the colour and angular momentum 
states of the QQ pair. Non-perturbative QCD effects, mediating the transition 
to the colour- singlet bound state H{J^^) containing the QQ pair are assumed 
to be first universal and secondly negligible for the dynamics of H (a/s, p^, etc. 
dependence). The normalization factors for the various states are not predictable, 
but once fixed phenomenologically, the model is (surprisingly) successful. 

The factorization assumption (|^) of the colour-singlet model, on the other 
hand, says that all non-perturbative effects are contained in a single term that 
can be expressed as the non-relativistic wave function of the bound state. In turn, 
relative production rates of different quarkonium states can be predicted. More- 
over, different states may have different dynamical dependences since only specific 
short- distance cross sections contribute to each state. However, the colour-singlet 
model fails in two respects. First, predictions for 5- wave states often are way off, 
and second logarithmic infrared divergences spoil the factorization in the case of 



P-waves, cf. section 2.1 



This failure of the colour-singlet model can be traced back to that of the un- 
derlying quark potential model. Relativistic corrections are essential for a descrip- 
tion that is both consistent for P-wave states and successful for S'-wave states. 
Recently, a rigorous QCD analysis of the annihilation decays of heavy quarko- 
nium has been presented based on recasting the analysis in terms of HQQET, an 
effective field theory designed precisely for this purpose [^. It allows the sep- 



aration of long and short distances, where the short- distance contribution may 
be evaluated perturbatively, i.e. as series in asijn). The long-distance part is 
parametrized in terms of matrix elements, which are organized into a hierarchy 
according to their scaling with 1/m. 

A similar analysis has been performed in the context of non-relativistic quan- 
tum chromodynamics (NRQCD) . In this approach, the calculations are orga- 
nized in powers of the average velocity of the heavy (anti-)quark in the meson 
rest frame. Contributions of different orders in v are separated according to the 
"velocity-scaling" rules. This factorization formalism has been extended to the 
production cross section of heavy quarkonium in processes involving momentum 
transfers of order m or larger. 

The inclusive production cross section of a quarkonium state H and the parton 
fragmentation function into H take the form: 

c=l,8d=6,8,... X '''' 

D^Uz) = E '^^'l^ia'^^ 0|^f(^>-^;A)|0). (15) 

c=l,8d=6,8,... X '''' 
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Equation (|T5|) expresses the cross section (and analogously the fragmentation 
function) as a sum of terms, each of which factors into a short-distance coeffi- 
cient F„(A) and a long-distance matrix element (0|C^(A)) (where n = {c,d,X}). 
The coefficients F„ are proportional to the rates for the production of on-shell 
heavy quarks and antiquarks from initial-state gluons and light quarks, and they 
can be computed as perturbation series in as{m). They depend on all the kine- 
matical variables of the production process. The matrix elements {0\O^\0) give 
the probability for the formation of the quarkonium state H from the QQ pair of 
state n, and can be evaluated non-perturbatively using, for example, QCD sum 
rules or lattice simulations. The dependence on the arbitrary factorization scale 
A in ([T5|) cancels between the coefficients and the operators. 



The expansion ([T^) is organized into an expansion in powers of v . Only a 
finite number of operators contribute to any given order in v"^. Since the coeffi- 
cients Fn are calculated as perturbation series in as{m), eq. ( [T3| ) is really a double 
expansion in agim) and v"^. For heavy quarkonia, the two expansion parameters 
are not independent: v ~ agimv) > as{m). Hence corrections of order v"' must 
not be neglected compared to those of order a;"(m). Which terms in (|15]) actually 
contribute to the production of a quarkonium \H) depends on both the as{m) 
expansion of Fn and the f ^ expansion of the matrix elements. The latter is de- 
termined by the variables that specify a given operator 0^{d,X). These are the 
dimension d of the operator, and two further variables c and X, which specify 
the quantum numbers of the heavy QQ pair in the Fock-state expansion of the 
heavy quarkonium H. 

The state of the QQ pair in a general Fock state can be labelled by the colour 
state of the pair, singlet (c = 1) or octet (c = 8), and the angular momentum state 
^^^^Lj of the pair denoted by X. The leading term of the Fock-state expansion is 
the pure QQ state (@) of the potential model. Higher Fock states are suppressed 



by powers of v, which follow from the velocity-scaling rules |Q. For example, 
if the QQ in the dominant Fock state \QQ) has angular- momentum quantum 
numbers ^'^"'"^Lj, then the Fock state \QQg) has an amplitude of order v only 
if the QQ pair has total spin 5* and orbital angular momentum L + 1 or L — 1 
(El transition). Necessarily, the QQ pair must be in a colour-octet state. The 
general Fock-state expansion therefore starts as 

|ff(nJ^^)) = 0{1)\QQC'^'Lj,1)) 

+ 0{v) \QQr^\L±l)j,,8)g) 

+ 0{v') \QQC'^'Lj,8) gg) + ... 

+.... (16) 

In the case of J/ production, for example, one finds up to and including the 
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order v"^ 




m 



(0|Or(6,35i;A)|0) 



(0|Or(8,35i;A)|0) + O(i;^) . 



(17) 



Upon dropping the O(f^) contribution and identifying (0|Oi (6,^S'i;A)|0) = 
3iVc|i?i5'(0)p/(27r), eq. (0) reduces to the famihar factorization formula (§) of 
the colour-singlet model. Generally, the standard factorization formulas of the 
colour-singlet model, which contain a single non-perturbative parameter, are re- 
covered in the case of S-waves at leading order in v"^. 

However, the factorization formula is the sum of two terms in the case of 
P-waves. In addition to the conventional term, which takes into account the 
production (or annihilation) of the QQ pair from a colour-singlet P-wave state, 
there is a second term that involves production (annihilation) from a colour-octet 
S-wave state. For example, the gluon fragmentation function into Xcj states 



The state \Xc.j) consists predominantly of the Fock state \QQ)., with the 
QQ pair in a colour- singlet ^Pj state. It also has an amplitude of order v for the 
Fock state \QQg), with the QQ pair in a colour-octet ^5*1, ^Di, ^D2, or ^D-^ state: 



\XQj) = 0{l)\QQ{'Pj, 1)) + Oiv)\QQCSu8)g) + 0{v)\QQCDj,,8)g) + Oiv') . 



The Fock state \QQ) contributes to the production at leading order in through 
the dimension-8 operator Oi"' {8,^Pj). The Fock state \QQg), with the QQ pair 
in a colour-octet ^Si state also contributes to the production at the same order 
in f^, through the dimension-6 operator Og"'' {6, ^Si), because the latter scales as 
relative to Of^(8,3Pj). 

Equation ([18|) , together with the evolution equations for the matrix elements 
contain the solution to the problem of the infrared divergence encountered in 
At leading order in as{m), the coefficient df\z;^Pj; X) in (pT^) (= d[''\z,X) in 
(||,|D) depends logarithmically on the factorization scale A, while df\z; ^Si; A) = 
ds{z) in (|,|l3) is independent of A. In these coefficients, A plays the role of an 
infrared cutoff. The A-dependence of the short- distance coefficients cancels the 
A-dependence of the long-distance matrix elements, for which A plays the role of 
an ultraviolet cutoff. 



(J = 0,1, 2) is|l,cf. (0,1): 





1 y'' {0\O^^-'(8,'Pj.X)\0) 

+ d^i\z;'S,;X){0\Or'{Q,'S,;X)\0) + 0{v') . (18) 



(19) 
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At leading order in v"^ and Q;s(m), the dimension-8 matrix element Oi(A) = 
(0|(9f"'(8, ^Pj; A)|0) is renormalization-scale-invariant, while the dimension-8 ma- 
trix element 08(A) = (0|(98 "'(6, '^5'i; A)|0) has a non-trivial scaling behaviour [Q: 



^ ^ /^^ 16 (2J+ l)Ol , , 



With the help of ( pOD and (§,0) it is straightforward to show that 

A^=^ = . (21) 
dX 

This is in accordance with the general expectation that physical quantities, such 
as fragmentation functions, are renormalization-group invariants, i.e. independent 
of the arbitrary factorization scale A. (In practice, yet, it might not always be 
easy to know at which value of A the matrix elements are evaluated.) 

To leading order in the evolution equations (|20|) can be solved analytically 
with the result 

27po as[m) 

This solution may be used to provide an order-of-magnitude estimate of the 
colout-octet matrix element Og in terms of the colour- singlet matrix element Oi 
by assuming that 08(A) can be neglected compared to the second term in (^) 
for as(A) = 1. Recall that, to leading order in f^, Oi = 9|i?p(0)p/27r. Note also 
that, solving (pO|) for constant a^, one recovers the approximate solution ([To|). 

The new factorization approach provides a consistent framework for the cal- 
culation of P-wave production and decays; Xcj decays can well be described and 
a reasonable value of ^^(m) is found |TP[. Detailed comparisons of Xcj production 
with experiments have not yet been performed. Based on (|18D or (P) one expects 
Xcj production rates to be approximately proportional to (2 J + 1). Preliminary 
CDF data |^ on the Xci to Xc2 ratio at high p^, where gluon fragmentation is 
the dominant production mechanism, indeed seem to confirm this expectation. 

The comparison is less fortunate for the total production rates measured at 
fixed-target energies. Experimentally, the ratio of Xci to Xc2 production is about 
2 : 3 ||19[. This still seems to hold at -^i sufficiently large for the gluon fu- 



sion to dominate over quark-initiated processes. However, one then expects, 
theoretically, <y{xci) / <^{Xc2) = in leading order in v"^ and ^^(m): The Landau- 
Yan theorem forbids Xci production through the fusion of two gluons, hence the 
colour-singlet matrix element (0|O^''^(8, '^Pi)|0) does not contribute at 0{a1). 
However, there is no contribution either from the colour-octet matrix element 
(0|C8''^(6, ^S'i)|0): two-gluon fusion into a coloured ^5*1 state is absent, since in 
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the non-relativistic limit the corresponding amphtude is equal to zero. Hence 
corrections to the ratio are truly of 0{as{m)) or of O(f^), i.e. there is no O(f^) 
correction that is enhanced by a short- distance factor l/as{m). (Such a correc- 
tion does, however, exist in the case of qq annihilation.) 

Similar problems exist in the description of J/ip (and tp', collectively denoted 
by ip in the following) production and decays. Up to and including the O(f^) its 
Fock-state expansion is 

mi-)) = 0{l)\QQCSur))+ 0{v)\QQCPj,8)g) 

+ Y: {0{v')\QQCS„c)gg) + 0{v')\QQ{'Dj,c)gg)} 

c=l,8 

+ Oiv')\QQCSo,8)g) + 0{v') . (23) 

Hence corrections to the predictions of the colour- singlet model are truly of 0{v'^), 
i.e. to 0(1) the new factorization approach coincides with the colour-singlet model 
that is known to be way-off the data. Large relativistic corrections are therefore 
needed. 

In the context of a meaningful expansion in powers of v"^, large relativistic 
corrections can only be accommodated if contributions, suppressed by power(s) 
of become important because the leading contribution is suppressed by addi- 
tional power (s) of as{m). However, the 0(i;^) corrections are not enhanced by a 
short-distance factor l/as{m) since they are genuine O(f^) corrections to the pro- 
duction of a colour-singlet ^Si state, cf. (|T^. Only at 0{v'^) does such an enhance- 
ment occur: a factor l/as{m) in gluon-gluon fusion to J/ijj via (0|Og (6, ^5*0)10) 
and (0|Og (8, ^Po,2)|0) (relevant for ip production at fixed-target energies), and a 
factor l/aKm) in gluon fragmentation into ip via (OlOg (6, ^S'i)|0). Treating the 
latter matrix element as a free parameter, the description of high-p^ ^' produc- 
tion at the Tevatron can indeed be rescued |]3^, but clearly more work is needed 
before the new factorization formalism is established as (the) successful theory 
of quarkonium production. 

Further work can come from three sources. First, further phenomenological 
studies are needed. Ideally, one would like to have a global analysis of the data on 
charmonium production from all high-energy processes. Comparisons with bot- 
tomonium production should confirm that relativistic corrections indeed decrease 
as the heavy-quark mass increases. 

Secondly, more theoretical analyses of the foundation of the factorization 
formalism of NRQCD in (full) QCD are needed. A key ingredient is the 
velocity-scaling rules, both for the ordering of operators of a given dimension and 
the Fock-state expansion. Using perturbation theory, each additional gluon asso- 
ciated with the (assumed) dominant QQ pair is ascribed an extra power of v via 
the identification v ~ as{mv), valid for a colour- Coulomb potential. Although 
this estimate may be underlined by the multipole expansion, we do not know of 
any rigorous derivation. One way would be the extension of the factorization of 
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HQQET from quarkonium decays to their production. Additionally, impli- 
cations of spin symmetry for the production (and decays) of heavy quarkonia 
should be investigated further. 

Thirdly, more experimental data will improve our understanding of bound 
states of heavy quarks. Concerning their production, one would in particular like 
to see ratios of the production rates of various charmonium states; as a function 
of pt at the Tevatron and HERA, and their xp dependence in hadro- and pho- 
toproduction. Interesting would also be the observation of a spin alignment of 
the heavy quarkonium states. Last but not least, valuable information will come 
from charmonium production in e^e~ collisions through improved measurements 
of 6-decays into charmonium, observation of a fragmentation contribution, and 
(at LEP2) charmonium production in two-photon collisions. 
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Figure 1: Lowest-order contributions to heavy-quark hadroproduction. 



Figure 2: Cross sections for b and c production in vrA^ collisions. For details see 
text (from ref. [||). 

Figure 3: Total charm production cross sections from pp and pA measurements 
compared to calculations. The curves are: MRS D— ' nic = 1.2 GeV, /i = 2mc 
(solid); MRS DO' = 1.2 GeV, /i = 2m^ (dashed); GRV HO = 1.3 GeV, 



= rric (dot-dashed); GRV HO rric = 1.5 GeV, fi = rric (dotted) (from ref. ||T0|). 

Figure 4: The 6-quark cross sections at CDF. For details see text (from ref. 0). 

Figure 5: Preliminary CDF data for prompt J/tp and tp' production compared 
with theoretical predictions of the total fragmentation contribution (solid curves) 
and the total leading-order contribution (dashed curves) (from ref. 



Figure 6: The ratio of (xci + Xc2) ^ J/ip to total J/ijj production as a function 
of cm. energy ^/s, by proton (open symbols) and pion beams (solid symbols) 
(from ref. |30|). 



Figure 7: The ratio of ip' to J/ip production as a function of cm. energy -^i, on 



proton (circles) and nuclear targets (squares) (from ref. ||30[| ). 



Figure 8: The ratio of ip' to J/ip production as a function of transverse momen- 



tum; the shaded strip shows the average value of Fig. 7 (from ref. ^0 



Figure 9: The differential J/ip production cross section {da[pN J/ipXj/dy) = 
2.5 X 10^^ {da[cc]/dy) aX y = 0, calculated with MRS D— ' PDF, compared with 
data (from ref. ||30|| ). 

Figure 10: Energy dependence of T production in pN collisions using the MRS 
D— ' PDF. Also shown (CR) is a phenomenological low-energy fit (from ref. ||30|| ). 

Figure 11: The J/ip longitudinal momentum distributions compared with pN 
(top) and pN (bottom) data using two parametrizations of the PDF, MRS D— ' 
(solid) and GRV (dashed) (from ref. g). 
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Figure 1: Lowest-order contributions to heavy-quark hadroproduction. 
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Figure 2: Cross sections for h and c production in TriV collisions. For details see 
text (from ref. [8]). 
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Figure 3: Total charm production cross sections from pp and pA measurements 
compared to calculations. The curves are: MRS D— ' rric = 1.2 GeV, fi = 2mc 
(solid); MRS DO' = 1.2 GeV, = 2me (dashed); GRV HO = 1.3 GeV, 
fi = rric (dot-dashed); GRV HO rric = 1.5 GeV, fi = rric (dotted) (from ref. [10]). 




Figure 4: The 6-quark cross sections at CDF. For details see text (from ref. [8]). 
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Vs = 1.8 TeV |t?^|<0.6 




Pt (GeV) 




Figure 5: Preliminary CDF data for prompt J/^ and ^' production compared 
with theoretical predictions of the total fragmentation contribution (solid curves) 
and the total leading-order contribution (dashed curves) (from ref. [26]). 
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Figure 6: The ratio of [xci + Xc2) — ^ J/i^ to total J/il^ production as a function of 
cm. energy y^, by proton (open symbols) and pion beams (solid symbols) (from 
ref. [30]). 
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Figure 7: The ratio of to J/il^ production as a function of cm. energy y^, on 
proton (circles) and nuclear targets (squares) (from ref. [30]). 
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Figure 8: The ratio of ^' to J/^ production as a function of transverse momen- 
tum; the shaded strip shows the average value of Fig. 7 (from ref. [30]). 
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Figure 9: The differential J/'^ production cross section [da[pN — )• J/il^X]/dy) = 
2.5 X 10-2 {da[cc]/dy) at y = 0, calculated with MRS D-' PDF, compared with 
data (from ref. [30]). 
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Figure 10: Energy dependence of T production in piV collisions using the MRS 
D— ' PDF. Also shown (CR) is a phenomenological low-energy fit (from ref. [30]). 
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Figure 11: The J/^ longitudinal momentum distributions compared with pN 
(top) and pN (bottom) data using two parametrizations of the PDF, MRS D— ' 
(solid) and GRV (dashed) (from ref. [30]). 
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